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Abstract 

PreUminary results from the BABAR experiment on charmless B decays to charged pions or 
kaons, and the measurement of the B^J/tpK* decay amplitudes are presented. The data sam- 
ple, collected at the asymmetric-energy i?-factory PEP-II at SLAC, comprises a total number of 
22.7 million T{AS) decays, corresponding to an integrated on-resonance luminosity of approxi- 
mately 21 fb"^ We measure the following CP-averaged branching fractions: B{B^^t:~^tt ) = 
(4.1 ± l.O(stat) ± 0.7(sys)) x 10-^ B{B^^K+tt-) = (16.7 ± 1.6(stat) l}:?(sys)) x 10-^ and an 
upper limit of B{B^^K^ K^) < 2.5 x 10^^, at 90% confidence limit. The measurement of the 
J/ipK* decay amplitudes results in R± = 0.160 it 0.032(stat) it 0.036(sys), and reveals a dominant 
longitudinal component. The phase of the longitudinal amplitude shows evidence for non- vanishing 
final-state interaction. 
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1 Charmless Hadronic B Decays 



Hadronic B decays to charmless final states provide a rich environment for CP violation studies as 
well as searches for new physics. Any b uud transitions, where the hadronic final state may be a 
CP eigenstate or not, represents a possible source for measuring the angle a = arg[— VtdV^^/KjrfFJJ,] 
of the unitarity triangle (UT). Indirect CP violation in B^B^ mixing results from the interference 
between direct and mixed decays. The weak phase a can be determined from the measurement 
of the time-dependent asymmetry of iJ^^vr+vr^, together with constraints from isospin conserva- 
tion Q. Moreover, a three-body Dalitz plot analysis of B^pn decays allows one to simultaneously 
fit for a and the contributing strong tree and penguin phases Direct CP violation arises from 
the interference between tree and penguin amplitudes contributing to the hadronic final state. It 
can be measured in charged as well as neutral B decays. Several charmless modes are self-tagging 
which facilitates the search for direct CP violation. Theories employing perturbative QCD Q or 
non-leading Factorization Approximation (FA) |Q] have been shown to provide predictions of charm- 
less branching fractions of tree and penguins dominated B decay modes. Their measurement can 
exhibit tests of the theory and may eventually be used to determine the UT angle 7 = arg(l^*j) ||5|. 
Finally, the B^^K~^ K~ branching fraction provides a measure of contributions from W exchange 
and hard rescattering topologies which are suppressed in the FA. 

Preliminary BABAR results on charmless B decays using approximately half of the present data 
set have been presented in Ref. Q. Significant signals have been found for the modes: B^^t:~^'k~ , 
B^^K^7r~ , B^^p^iT~ and B^^-q'K. Updates of these analyses as well as measurements of a large 
number of new charmless modes are being prepared within the BABAR Collaboration. Important 
activity is also invested into the analyses of 6^57 decays for which results are forthcoming. 

In this note we will concentrate on the measurement of the charged two-body modes B^^h^h'~ , 
with /i^') being pion or kaon. We use (almost) the full year 1999-2000 data set which amounts to 
(22.57 lb 0.36) X 10^ BB pairs, corresponding to an integrated luminosity of approximately 21 fb~^ 
collected by the BABAR detector at the PEP-II storage ring. Background studies utilize approx- 
imately 3 fb^^ of data taken below the BB threshold. A "blind" analysis methodology has been 
adopted throughout, so that the signal region for each mode remained hidden until finalization of 
event selection and systematic studies. 

1.1 B^^h+h'- Candidate Selection 

We impose efficient track quality requirements on charged tracks which must have a minimum 
transverse momentum of lOOMeV/c in the laboratory frame (LAB). Charged pions and kaons are 
identified by measuring the Cerenkov angle, 9c, of photons produced while traversing a radiative 
medium made of synthetic fused silica, the DIRC |^ , surrounding the drift chamber in the barrel Q . 

We select B^^h'^h'~ candidates based on the energy-substituted mass, tit-es, where ^/s/2 {^/s being 
the center-of-mass energy) is substituted for the candidate's energy, and the difference IS.E between 
the 5-candidate energy and \fsl1- Preselected B candidates are required to satisfy 5.2 < m-Es < 
5.3 GeV/c^ and |A£^| < 0.42 GeV. The dominant background for all charmless modes stems from 
continuum qq production, which exhibits a jet-like structure that distinguishes it from the more 
spherical topology of generic BB events. To suppress this background we use a Fisher discriminant, 
J- , combining the charged and neutral momentum flow from the rest of the event relative to the 
candidate's thrust axis, characterized by nine concentric, 10° wide cones Q. 
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Figure 1: Energy substituted mass spectra. Left: B tt signal with single Gaussian fit. Right: 

background from AE sidebands of B^h^h'^ candidates, fit to an ARGUS shape function [17|. 



1.2 Calibration of the Discriminating Variables 

The signal and background distributions of the discriminating variables are controlled by inde- 
pendent calibration samples obtained from: /S.E and mss sidebands and off-resonance data to 
study background characteristics, exclusive B~ ^D^'k~ ^{K~ 'k^)t:~^ decays to gauge the AE 
and mES resolution of the Monte Carlo (MC) simulation with respect to the data, inclusive 
D*^ ^D^t:^ ^{K^ it^)t:^ decays to measure the particle identification (PID) performance of the 
DIRC, and finally MC simulated signal and background events. The calibration measures a rriES 
resolution of 2.6 MeV/c^ for signal events (see Fig. [l| for signal and background distributions), 
and a AE resolution of approximately 26 MeV. The latter value has been rescaled to account for 
differences between data and MC simulation with a conservative error estimated to be ±5 MeV. 
The output of the Fisher discriminant is found to be well reproduced by the MC simulation. A 
sample of I?*-tagged D^^K~t:~^ decays is used to parameterize the 9c distributions for pion and 
kaon tracks as a function of the track's polar angle. The K/tt separation varies from 2 to more 
than 8 standard deviations across the relevant momentum range between 1.7 and 4.2 GeV/c in the 
LAB. 



1.3 Likelihood Fit 

Signal yields in all three h^h'^ modes are determined simultaneously from an 8 parameter, unbinned 
maximum likelihood (UML) fit, incorporating niES, AE, T , and the measured Oc for each track. 
All candidates passing the preselection are included in the fit. The probability density function 
(PDF) for an event z, with a number of signal (background) candidates, A''(b)/ife', is given by 

Pi(/x)oc (A^w7'f'(/^)+ iVbwn''''''' (/"))+ asymmetries (1) 

h<h'=-K,K 



with (fi) = {mES: ^E,T,9'^ ,6^, ) being the discriminating variables that enter the UML fit. The 
combined PDF, V'^^^^^' {jj) , is the product of the PDF's of the individual variables. The total 
extended log-likelihood [|lO| for candidates, which is to be maximized in the fit, reads 

N 

C = Y.V.{^l) - N' , (2) 

1=1 

where A^' = J2h h'i-^h,h' + ^hhh') is the Poissonian extension of the likelihood. A bias, enhancing 
the contributions from 7r"'"7r~ and K^K~ with respect to K~^Tr~, has been identified in the 6 
parameter UML fit previously applied ||6|. This was due to the implicit assumption of randomly 
paired background tracks. The addition of 2 independent background parameters provides the 
remedy to the problem. 



1.4 Systematic Uncertainties 

Systematics uncertainties arise from the imperfect knowledge of the PDF shape parameterizations. 
Associated errors are quantified by varying the parameters within their statistical errors and to 
cover disagreements between MC and data. Alternate parameterizations are used in addition. A 
large number of cross checks have been performed to validate the fit procedure. Toy MC studies 
showed no bias in the fit outputs. A cut-based analysis, performed in parallel, agrees with the 
results from the UML fit. 



1.5 Results 

With signal efficiencies of e^^ = 0.448±0.008, exn = 0.447±0.008 and crk = 0.429±0.008, we find 
signal yields A^^^ = 41 ± 10 ± 7, Nk-jt = 169 ± 17lJ^ and Nkk = 8.2tli ± 3.3, from the likelihood 
fit, where the uncertainties are statistical and systematic. The correlations between these numbers 
are below 15%. The systematic errors are dominated by the uncertainties in the description of J- 
and the resolution on AE. Figure ^ shows the uies distributions for candidates of the three final 
states from the cut based analysis. 

The branching fractions are obtained from the signal yield via B{B^h^h'^) = Nfih' /i^kh'^Bs')^ 
giving 

B{B°^TT+TT~) = (4.1 ± 1.0 ± 0.7) X 10-6 , (3) 
B{B^^K^TT~) = (16.7 ± 1.6l{:?) X 10"^ > (4) 
BiB^^K+K-) < 2.5 X IQ-^ (90% CL) , (5) 

where the uncertainties are statistical and systematic. The statistical significance is 4.7 (15.8) 
standard deviations for the 7r~^TT~ {K~^7r~) final state. We do not find a significant signal yield 
in the mode B^^K^ K~ and hence quote an upper limit at 90% confidence level. The results 
are consistent with our previous preliminary numbers B , based on less than half the present data 



set. A comparison with published branching ratios from CLEO |11] and preliminary results from 



BELLE [12 1 is given in Fig. y, showing consistency. 
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Figure 2: Distribution of rriEs for the final states vr+vr (upper plot), K+vr (middle plot) and 
K^K~ (lower plot), from the cut based analysis. 



2 Measurement of the B^J/i/jK* Decay Amplitudes 



The decay B^^J/t/j K*^ , K*^^KgiT^ can be used in time dependent CP asymmetry studies to 
measure the CP violating parameter sin2/3. In contrast to the pure CP-odd eigenstate B^^J/ip Kg, 
this pseudoscalar to vector-vector decay gives a mixture of CP-odd and CP-even eigenstates since 
orbital angular momenta L = 0, 1,2 of the {J/iIj ,K*^) system are allowed. If one CP contribution 
dominates, sin2/3 can be measured as in the B^^J/ipK^ case, using the time distribution only. 
Otherwise, both the angular and time distributions are needed to avoid a partial cancellation of 



the observed CP asymmetry |13]. Assuming the weak decay to be isospin invariant, the amplitudes 
can be measured using all K*^KTr modes, which greatly increases statistics. 

Three amplitudes are necessary to describe the decay. In the transversity basis |l4|, used in 

the amplitudes Au, Aq and A±, have definite CP 



previous analyses by CLEO |15] and CDF ||Tf^ 
eigenvalues: the latter is CP-odd and the former two are CP-even. The longitudinal polarization 
fraction of the decay is measured by \Aq\'^ = F^/r, whereas and A± govern the transverse 
polarization and correspond to parallel and perpendicular polarizations of the two vector mesons, 
respectively. Hence, measures the P wave contribution to the decay. The measurement 

of the decay amplitudes also provides a test of the factorization hypothesis. In the factorization 
scheme, the weak decay is described by a product of two hadronic currents, J/ip and B^K*, where 
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Figure 3: Comparison of measured branciiing ratios to h^h' final states. The shaded bands give 
the averages for the significant measurements. 



final-state interaction is neglected. If factorization holds, the corresponding amplitudes should have 
trivial relative phases of or vr. 

We present preliminary results on a complete angular analysis in the transversity basis considering 
the decays B^^J/tpK*^ and B^^J/ipK*~^, where the K*^ and K*~^ are reconstructed in the modes 
KgTT^, K~^7r~ and Kg7r~^ and K~^Tr^, the Kg decaying into two charged pions, and the J/ijj through 
decays into e~^e~ and fJ-^fJ,^. The data sample used corresponds to the full year 1999-2000 recording 
of (22.7 ± 0.4) X 10^ BB pairs (see Section [l|). 

2.1 B^J/^K* Candidate Selection 

Charged tracks are used in a fiducial polar angle range where PID efficiencies are well understood. 
Electrons are identified requiring their momentum to be compatible with their energy deposition 
in the electromagnetic calorimeter (EMC), and their dE/dx to be compatible with the electron 
hypothesis. Muons are required to penetrate at least 2 interaction lengths in the detector and to 
leave a low number of hits per layer in the Instrumented Flux Return (IFR). If a muon penetrates 
the EMC, it is required to deposit energy consistent with a minimum ionizing particle. Kaon 
candidates are required to pass a pion veto based on DIRC and dE/dx information. 

J/ip candidates are selected requiring both leptons to be identified. The vertexed lepton pair is 
required to have an invariant mass in the range 3.06-3.14 GeV/c^ for muons and 2.95-3.14 GeV/c^ 
for electrons. Kg are reconstructed as vertexed pairs of charged tracks with an invariant mass 
between 489 and 507 MeV/c^ and a flight direction which is compatible to originate from the 
interaction point. Photons are deflned as neutral clusters with an energy greater than 30 MeV and 
photon-like shower shapes. Neutral pions are selected from pairs of photons with an invariant mass 
between 106 and 153 MeV/c^. The J/V', K° and vr'' are fitted to their respective nominal masses. 
K* are formed as Kn combinations with an invariant mass within ±100 MeV/c^ of the nominal 
K*{892) mass. 



Finally, B mesons are formed by combining J/ip and K* candidates. Furthermore, requirements 
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Figure 4: Background-subtracted i^vr mass distribution in the mode B^^J/t/^K*^, K 



*o 



on the cosine of the K* hehcity angle reduce the cross feed (CF) from J/tl) (K-jt^)* modes, where 
the vr^ is lost, and the self cross feed (SCF) due to a wrongly reconstructed vr'' picked up by 
the reconstruction. The (S)CF is the most important background and accumulates in the tti-es 
signal region. B meson candidates are selected with AE (required to be between —70 MeV and 
+50 MeV for channels involving a vr'^ and within ±30 MeV in the other cases) and tit-es- For the rriEs 
signal region between 5.27 and 5.29 GeV/c^ we obtain reconstruction efficiencies of 9.9%, 23.9%, 
17.2% and 13.8% for the J/ipK^Tr^, J/tpK+ir-, J/tljK^iT+ and J/V^K+7r° modes, respectively The 
corresponding (S)CF contamination amounts to (1.6)15.8%, (0.6)2.4%, (0.5)3.0% and (2.2)15.7% 
of the signal contributions. The signal yields are obtained by a fit of the mss spectra describing 
the background by an ARGUS function |17] and the signal by a Gaussian: [J/i/j K^tt^) = 39 ± 6, 
N (J/ipK+TT-) = 530 ± 32, N {J/ipK^7r+) = 131 ± 12 and N {J/ip K+tt'^) = 195 ± 15. 



As shown in Ref. |15], Ktt production other than via the i^r*(892) contributes to B^J/'ijj Kt: final 
states. Figure Q shows the Ktt mass distribution for the B^J/if) K^tt^ mode, after subtraction 
of combinatoric Kt: background from the low mES region. We observe clear evidence of resonant 
excitation around 1430 MeV/c^. 



2.2 Angular Distribution and Fitting Method 

The transversity frame is defined in the J/ip rest frame. The K* direction defines the —x axis. 
The Ktt plane defines the (x,y) plane, with y such that Py{K) > 0. The z axis is normal to that 
plane so that the {x,y,z) frame is right handed. The transversity angles ^tr and (ptr are defined 
as the polar and azimuthal angles of the originating from the J/V' • The K* helicity angle 9k* 
is the angle between the K direction and the direction opposite to the J/V' in the K* rest frame. 
Defining: 

fi = + 9/(327r) • 2 cos^ Ok* (1 - sin^ 0^^ cos^ </>tr) 
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the angular distribution g{cos9ij.,cos9K* 



reads 
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d^r 



9 = 



r dcOS ^tr dcOS 9k* d(/)tr 

2 , f I /l..|2 , ^ I /I , |2 



/il^or + M^iir + /3|^±r + /4lm(A||*A^) + /5Re(Ao*A||) + /6lm(Ao*Aj 



(6) 



The parameter is +1 for and B'^' and —1 for i?^ and S*^. For the CP eigenstate K^tt^, a B^B^ 
from the decay of an T{AS) has C = ±1/(1 + 2;^). However, the initial B^ flavor is not determined 
in this analysis so the appropriate value for C is the average, namely zero. 

The amplitudes are fitted by means of a UML fit exploiting the angular and tit-es information. All 
B candidates within 5.2 < mES < 5.3 GeV/c^ enter the fit, where the ttt-es sidebands are used to 
normalize non B^J/ipK* background. The biases on the amplitudes due to limited acceptance 
and background are corrected considering only the projections on the basis of the /j's |18]. 

Signal events are described by the PDF 

5obs = gi^^j) X e{uj)/{e) , (7) 

where tOj stands for the three angles of the observed event j, and e{ujj) is the efficiency for this event. 
Rewriting Eq. (|) as g = ELi where the = 1,6) are |^oP, Im(A||*A^), 
Re(Ao*A||) and Im(Ao*A_L), the mean efficiency (e) reads 



6 

2 = 1 



(8) 



where the ^.i = J fi x e are amplitude independent. The signal part of the log- likelihood 

Nobs 

>Csignal = log (fi'obs('2j)) , 

where A^obs is the number of observed events, becomes 

Nobs Nobs / 6 \ 

^signal = (ai^i)) + J2 log (^('^i)) - ^obs log ( ^ A^^i ) • 



(9) 



(10) 



\i=l 



Since the e{u}j) do not depend on the amplitudes, the term J2f=i' log(e('^j)) is left out: detailed 
knowledge of the acceptance is not necessary in this approach. 

The coefficients are evaluated using a detailed simulation of the BABAR detector. Separate series 
of ^i's are used for each B^J/jpK* channel. The values found for ^j, {i = 1,3), are close to the 
efficiency of the final state, with a lower value for ^1, especially in modes involving tt^, due to the 



Table 1: Systematic uncertainties. Details are given in the text. 













<t>\\ 


MC Stat. 


0.006 


0.006 


0.007 


0.04 


0.06 


Backgr. 


0.002 


0.005 


0.006 


0.01 


0.00 


Track. & PID 


0.002 


0.006 


0.004 


0.00 


0.02 


Heavy K* 


0.005 


0.035 


0.031 


0.04 


0.02 


Total 


0.008 


0.036 


0.033 


0.06 


0.07 



Table 2: Preliminary amplitude moduli-squared and phases. The uncertainties are statistical and 
systematic. 
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requirement on cos 9k*- The ^i, {i = 4,6), are compatible with zero. 

The PDF of the non B^J/ip K* background g'B^{(jJj) is written as a distribution gs of amplitudes 
B^{i = 1,6): g'B^iujj) = gsi'^j) X e('^j)/(e)- The shape and level of the (S)CF background is 
amplitude dependent. The bias due to (S)CF background can be corrected by a modification of 
the ^j, obtained by integrating, in addition to the signal as before, over the (S)CF contributions. 
In contrast to the depend on the amplitudes used in the simulation. The maximum bias on 
the fitted amplitudes was found to be a few parts per thousand. 

Finally, the complete log-likelihood has the form 




where G(mEs) and B{mEs) are the Gaussian signal function and ARGUS background function, 
respectively. The parameter x mcasTircs the signal fraction. The normalization of the amplitudes 
is enforced by the likelihood extension J\f = A'obsd^oP + + |^±P)- 

2.3 Systematic Uncertainties 

The limited MC statistics (32000 events per channel) gives rise to a systematic uncertainty on 
acceptance and cross feed corrections. A fraction of the feed-across background may be absorbed 
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Figure 5: One-dimensional projections of the three transversity angles. The result of the fit is 
superimposed. The data are corrected for bin-by-bin acceptance. Top plot: modes without Tr*^ in 
the final state, and bottom plot: modes with -k^ . 



in the fit by the non B^J/ipK* background component. The size of this effect is evaluated 
on MC samples, where a set of fits are performed with the non-J/^j^K* background (ARGUS) 
function removed. The largest deviation found is taken as an estimate of the associated systematic 
uncertainty. Discrepancies between the true detector and the simulation affect the acceptance 
correction. The scale of the effect is estimated by varying the tracking and particle identification 
efficiencies at MC level. Non-resonant Ktt combinations, higher mass Ktt resonances and their 
interference with the X*(892) can contribute to the events in the K* mass window. The size of 
the potential bias from these events is estimated by the magnitude of the shift in the fitted angular 
amplitudes when events in the region \mK'w — 1-2| < 0.1 GeV/c^ are included in the fit. 

Table |l] quantifies the systematic effects of the contributions listed above on the fit observables. 
2.4 Results 

The UML fit results for amplitude moduli-squared and phases are given in Tab. H. The errors 
account for statistical and systematic uncertainties. One-dimensional projections of the fitted 
amplitude spectrum on the three transversity angles are shown in Fig. ^ for modes without vr^ 
(upper plots) and with vr^ (lower plots) in the final states. Figure ^ visualizes the quantitative 
effects of the different corrections on the amplitude moduli-squared. Only the acceptance correction 
exceeds the range defined by the statistical uncertainties. Figure |^ shows the comparison between 
this preliminary measurement and results from BELLE [l^ (preliminary), CDF |16| and CLEO [15|, 
which are found to be consistent. The shaded bands indicate the weighted averages. The BABAR 
values improve the precision by about a factor of two. 



We find a dominant longitudinal component and a small P wave contribution. If sin2/3 is measured 
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Figure 6: Results of the amplitude fits and effects of the corrections applied. The results are shown 
in the allowed domain of the amplitude moduli-squared defined by the conditions |^oP + + 
|A_Lp = 1 and \Ai\'^ > 0, (i = 0, ||,_L). Left: the entire allowed domain, right: zoom by a factor 
of ten into the (b, c, d) region, (a), angular fit with no acceptance correction or ttt-es dependence, 
(b), acceptance correction applied, (c), tues dependence added, and (d), (S)CF correction added. 
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Figure 7: Comparison of measured amplitude moduli-squared and phases (see text for references). 
The Belle fit assumed no final-state interaction. The shaded bands give the averages. 



from a time dependent CP asymmetry analysis in the mode B^J/i/j K^tt^ , the dilution introduced 
amounts to D± = (1 — 2|A_|_p) = 0.68 it 0.10. The uncertainty on contributes with relative 

14% to the total error of sin2/3. This measurement of implies a deviation from vr by 3 standard 
deviations (3.4 for the world average), and thus represents evidence for non- vanishing final-state 
interaction. 
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